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Microcanonical entropy inflection points: Key to systematic understanding 

of transitions in finite systems 
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We introduce a systematic classification method for the analogs of phase transitions in finite sys- 
tems. This completely general analysis, which is applicable to any physical system and extends 
towards the thermodynamic limit, is based on the microcanonical entropy and its energetic deriva- 
tive, the inverse caloric temperature. Inflection points of this quantity signal cooperative activity and 
thus serve as distinct indicators of transitions. We demonstrate the power of this method through 
application to the long-standing problem of liquid-solid transitions in elastic, flexible homopolymers. 
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Structure formation processes are typically accompa- 
nied by nucleation transitions, where crystalline shapes 
form out of a liquid or vapor phase. Thus, nucleation 
is governed by finite-size and surface effects. For small 
physical systems, it is difficult to understand thermody- 
namic transitions of this type, as they strongly depend 
on system size. 

Cooperativity refers to collective changes in a sta- 
tistically significant fraction of the degrees of freedom 
in a system, which transforms the system into a new 
macrostate. In the thermodynamic limit of an infinitely 
large system, the ensemble of macrostates sharing simi- 
lar thermodynamic properties would be called a "phase" 
and the transformation a "phase transition" . The de- 
scription of such a transformation in a finite system 
is more subtle, as it cannot be described in the tradi- 
tional Ehrcnfcst scheme of singularities in response quan- 
tities. However, statistical physics and thus thermody- 
namics are also valid for systems with no thermody- 
namic limit. Examples include the structure formation 
in small atomic clusters and all biomolccules. This is 
particularly striking for proteins, i.e., heterogeneous lin- 
ear chains of amino acids. The fact that the individ- 
ual biological function is connected with the geometrical 
shape of the molecule makes it necessary to discriminate 
unfolded (non-functional) and folded (functional) states. 
Although these systems are finite, they undergo a struc- 
tural transition by passing a single (or more) free-energy 
barrier(s). Since these finite-system transitions exhibit 
strong similarities compared to phase transitions, we ex- 
tend the terminology once defined in the thermodynamic 
limit to all systems exhibiting cooperative behavior. 

In this paper, we introduce a commonly applicable 
and simple method for the identification and classifica- 
tion of cooperative behavior in systems of arbitrary size 
by means of microcanonical thermodynamics |l| . It also 
includes the precise and straightforward analysis of the 
finite-size effects, which are important to a general un- 
derstanding of the onset of phase transitions. This is 
in contrast to canonical approaches, where detailed in- 



formation is lost by averaging out thermal fluctuations. 
Re-gaining information about finite-size effects in canon- 
ical schemes, e.g., by the investigation of the distribution 
of Lee- Yang zeros in the complex temperature plane [2j 
or by inverse Laplace transform [3|], is complicated. 

The identification of transitions is associated with a 
distinct definition of transition points such as a tran- 
sition temperature. In the canonical representation of 
finite systems these usually differ, e.g., peak structures 
of thermodynamic quantities, such as the specific heat 
and fluctuations of order parameters as functions of the 
heat-bath temperature. This makes it impossible to fix 
a unique transition point. In the microcanonical analy- 
sis, the temperature is defined via the curvature of the 
caloric entropy curve and thus all transition signals in the 
microcanonical entropy can be directly associated with a 
transition temperature. 

After introducing the method, we apply it to liquid- 
solid and solid-solid transitions occurring for elastic, flex- 
ible polymers, which have been under debate for quite 
some time. In contrast to the rather well-understood 
coil-globule collapse transition, the formation of highly 
compact crystalline, amorphous, or glasslike structures 
intricately depends on the precise relation of intrinsic en- 
ergy and length scales in the system |j49[. 

In recent work, the microcanonical analysis has suc- 
cessfully been applied in aggregation studies of coarse- 
grained polymer and peptide models, where a nucleation 
process was found to be an energetically ordered hier- 
archy of individual structural subphase transitions [lOJ ] . 
Caloric approaches have also been used to investigate the 
folding behavior of proteins [Ill4l5j and the structural 
phases of polymers with stiff bonds [8( , as well as poly- 
mer adsorption transitions 16l4l8|. Other applications 
include the formation of galaxies [19( , the cluste ring and 
fragmentation of atomic clusters and nuclei |llj20l 1 2 II , 
and order-disorder transitions in spin systems |l|, |3, 122 - 
24j . Most of these studies are aimed at using the micro- 



canonical analysis as an alternative approach to inves- 
tigating finite-size scaling properties. However, a sys- 



tematic scheme for the classification of transitions in 
the respective finite systems has remained lacking. The 
method introduced here closes this gap through the in- 
troduction of an Ehrenfest-like analysis based on micro- 
canonical entropy inflection points. 

A fundamental property of each physical system, and 
the central quantity for our method, is the microcanon- 
ical entropy S(E) = ks \ng(E), where g(E) is the den- 
sity of states for a given energy E (in the following, 
we will set ks = 1)- Alternatively, a volume entropy 
can be defined via the integrated density of states by 
S'(E) = k B \nG{E) with G{E) = J^dE'g(E') [HJ], 
which is virtually identical with g(E)AE (\AE/E\ < f ) 
in the transition regions |l0| . It has been argued that 
only G(E) is consistent with the classical equiparti- 



tion theorem [261 . [2 
much less obvious 




however, its physical meaning is 
Therefore, we will continue us- 
ing S(E) instead. It should also be mentioned that g(E) 
is the "natural" output provided, e.g., by generalized- 
ensemble Monte Carlo methods. Among the most promi- 
nent of these methods are multicanonical [28[ and Wang- 
Landau |29( sampling, which enable a precise numerical 
estimation of this quantity over hundreds or even thou- 
sands of orders of magnitude |g, [Z] . 

A qualitative change in the interplay of entropy and en- 
ergy in the system is signaled by noticeable alterations in 
the curvature of S(E), which are quantitative measures 
for the strength of cooperativity of the associated transi- 
tions. For finite systems exhibiting transitions withphase 
separation, S(E) can even possess convex regions [l|, & 1~ 
though it is a strictly concave function in the thermody- 
namic limit. In this case, the slope of a tangent at each 
point of the curve is unique, and it is common to define 
the reciprocal microcanonical temperature via the caloric 
derivative of S, 



{3(E)=T- 1 (E) = (dS/dE) Ny , 



(1) 



where system size N and volume V arc kept constant. 
In the thermodynamic limit, where fluctuations about 
the mean energy become negligible, the canonical and 
microcanonical ensembles are identical, and the canoni- 
cal (or heat-bath) temperature equals the microcanonical 
temperature. This is not the case for a finite system ex- 
periencing a structural transition, where different quan- 
tities vary in their fluctuation properties, rendering an 
identification of transition points impossible. Since the 
complete phase behavior is already encoded in S(E), it is 
useful to consider (3(E) as a unique parameter to identify 
transition points. 

We further propose to analyze the monotonic behav- 
ior of (3(E), expressed by its derivative with respect to 
energy, 



j(E) = d(3(E)/dE = d 2 S/dE 2 



(2) 
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This will allow for the introduction of a systematic clas- 
sification scheme of transitions in finite systems. In prin- 



FIG. 1: (Color online) Inverse temperature /3(e) and its 
derivative 7(e) as functions of the energy per particle, e = 
E/N , exemplified for an elastic polymer with 102 monomers. 
The maxima of 7(e) indicate transitions between the struc- 
tural phases A and B at e tr B and B and C at e B . The as- 
sociated points /3(e£ B ) = (3(f and /3(e t BC ) = /3 t BC define the 
transition temperatures T t f B = (/3£ B ) _1 and T t BC = (/3 BC ) _1 - 
According to our classification scheme, the transition between 
A and B is of second order, since the slope of the inflection 
point is negative. On the other hand, Bf>C is a first-order 
transition as the respective slope at /3(e B ) is positive. The 
non-monotonicity of /3(e) in this region, called "backbending" , 
is a typical signal of phase coexistence. The latent heat Aq 
is defined as the energetic width of this transition region. 

ciple, this can also be used for scaling analyses towards 
the thermodynamic limit. 

We define a transition between phases to be of first 
order if the slope of the corresponding inflection point of 
(3(E) at E = E tT is positive, i.e., 7 tr = 7(£tr) > 0. Only 
in this case is the temperature curve non-monotonic and 
there is no unique mapping between (3 and E. Physically, 
both phases coexist in the transition region. The overall 
energetic width of the undercooling, backbending, and 
overheating regions, obtained from a Maxwell construc- 
tion, is thus identical to the latent heat. Therefore, for 
a first-order transition, Aq > 0. In the case that the 
inflection point has a negative slope, 7 tr = j{E tT ) < 0, 
the phases cannot coexist and the latent heat is zero. 
In complete analogy to phase transitions in the thermo- 
dynamic limit, we classify such transitions as of second 
order. Since the inflection points of (3(E) correspond to 
maxima in j(E), it is therefore sufficient to analyze the 
peak structure of j(E) in order to identify the transition 
energies and temperatures. The sign of the peak values 
classifies the transition. This very simple and general 
classification scheme applies to all physical systems. 

Figure [1] illustrates the procedure for the identifica- 
tion of the transitions by means of inflection-point anal- 
ysis, where the inverse temperature (3 and its energetic 
derivative 7 are plotted as functions of the reduced en- 
ergy e = E/N, with N being the system size. As a 




FIG. 2: (Color online) Caloric temperature curves T(e) = 
/3 _ (e) for a selection of elastic, flexible polymers with chain 
lengths in the interval N = 13, ... , 309 (from right to left). 
Curves for chains with magic length (N — 13, 55, 147, 309) 
are bold. The relevant inflection points, indicating the con- 
formational transitions on the basis of our analysis method, 
are marked by • symbols. 

first example, we consider an elastic flexible homopoly- 
mer with N = 102 monomers. This system exhibits four 
structural phases [7J: two solid icosahedral phases (A: 
Mackay, B: anti-Mackay), a globular liquid phase (C), 
and the random-coil phase (D). In Fig. [TJ the transi- 
tions can indeed be uniquely identified (since the Cf>D 
transition occurs at much higher energy and tempera- 
ture, it is not included, but can also easily be found by 
inflection-point analysis; it is a second-order transition at 



= CD 



-1.21, P\f « 1.08). A l st -order liquid-solid tran- 
sition BoC is characterized by 7 BC > at ef r c « —4.35 
(/3 BC sa 2.97). The width of the energetic transition re- 
gion corresponds to the latent heat, Aq BC , which is ob- 
viously nonzero because of the backbending effect or the 
coexistence of both phases in this region. The 2 nd -order 
transition AoB is found at e£ B w -4.58 (/3^ B w 3.34) 
by an inflection point with negative slope (7^ B < 0). 

In order to demonstrate the capability of our method 
to systematically analyze all transitions in finite sys- 
tems, we estimate the transition points for the entire 
set of elastic Lennard- Jones homopolymcrs with N = 
13, ... , 309 monomers. In the liquid and solid regimes, 
the structural behavior of these polymers is very sim- 
ilar to rare-gas systems consisting of N atoms, which 
also form compact, crystalline clusters at very low tem- 
peratures 0, Hj, |3fJ. We employ the standard model 
for flexible, elastic polymers, where the monomers in- 
teract via a truncated-shifted Lennard- Jones potential, 
£lTV*j) = E hJ (min( rij ,r c ))-E h j(r c ) with E L ,(r tJ ) = 
4e[(<r/rij) 12 — (c/rij) 6 ], where r^ is the distance between 
two monomers located at r^ and r, (i,j = 1, ...,iV), 
and e = 1 and a = 2~ 1 ' 6 ro, with the potential min- 
imum at ro = 0.7 and the cutoff at r c = 2.5er. Ad- 
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FIG. 3: (Color online) Transition temperatures T tT (N) of con- 
formational transitions for small elastic polymers with chain 
lengths N = 13, ... , 309 in the liquid-solid and solid-solid 
transition regimes, obtained from inflection-point analysis. 
First-order transition points are marked by red symbols (•), 
second-order transition points by blue symbols (x). Also 
shown is a fit for the liquid-solid transition temperature to- 
wards the thermodynamic limit N — > oo (dashed line). 

jacent monomers are connected by finitely extensible 
nonlinear clastic (FENE) anharmonic bonds [31J, |32|, 
E FE NE(n l+1 ) = -KR 2 \n{l-[{r ll+1 -r Q )/RYY'\ The 
FENE potential minimum is located at r and diverges 
for r — > r ± R (in our simulations R = 0.3). The 
spring constant K is set to 40. The total energy of 
a polymer conformation X = (ri,...,rjv) is given by 

£(X) = Eli Ef= i+1 E^Hnj) + Eti 1 ^FEN E (r M+1 ). 

Figure [2] shows the caloric temperature curves for elas- 
tic polymers with various chain lengths in the liquid and 
solid regimes, calculated from highly accurate density of 
states estimates obtained in sophisticated multicanoni- 
cal Monte Carlo simulations [33j. The identified inflec- 
tion points associated with conformational transitions are 
indicated by • symbols. As expected, there is no gen- 
eral and obvious relation of the behavior of chains with 
slightly different lengths. This is due to the still dom- 
inant finite-size effects of the polymer trying to reduce 
their individual surfacc-to-volume ratio, which therefore 
strongly depends on optimal monomer packings in the in- 
terior and on the surface of the conformations. For exam- 
ple, for chains of moderate lengths (N < 147 [7|, |9(), the 
different behavior can be traced back to the monomer ar- 
rangements on the facets of icosahedral structures, known 
as Mackay and anti-Mackay overlayers [34[. Solid-solid 
transitions between Mackay and anti-Mackay structures 
are alsopossible under certain conditions in these sys- 
tems [7],|9|. This can be seen in Fig.[3l where all transition 
temperatures T tr (N) = /3 t ~ 1 (A r ) for liquid-solid and solid- 
solid transitions are plotted in dependence of the chain 
length N [35|. Symbols • indicate first-order transitions, 
which for N > 38 can be associated to the respective 



liquid-solid transitions, whereas symbols x mark second- 
order transitions. 

If the associated transition temperatures are smaller 
than the liquid-solid transition temperatures, the sym- 
bols indicating second-order behavior belong to solid- 
solid transitions, e.g., transitions between geometrical 
shapes with Mackay or anti-Mackay ovcrlaycrs. Note the 
different behavior for "magic" chain lengths iV mag i c = 
13, 55, 147, 309, . . ., in which icosahedral Mackay ground 
states form. Figure |3] also gives evidence for the conver- 
gence of the solid-solid and liquid-solid transition temper- 
atures when N approaches a magic length. This behavior 
repeats for each N interval that finally ends at a certain 
magic length iV mag i C , where both transitions merge into 
a single first-order liquid-solid transition. The influence 
of the solid-solid effects weakens with increasing system 
size, while the liquid-solid transition remains a true phase 
transition in the thermodynamic limit. Inserted into the 
plot is a fit function T tI (N) = T ls - aN" 1 ^ which sug- 
gests an estimate for the thermodynamic phase transition 
temperature T^ sa 0.64. 

Summarizing, we have introduced a general method 
for the analysis of phase transitions in small systems 
based on the central quantity of any statistical sys- 
tem, the microcanonical entropy, and applied it to the 
long-standing problem of structural transitions of flexible 
polymers. Advanced Monte Carlo simulation techniques 
such as multicanonical sampling |28[ and the Wang- 
Landau method [23 ] enable precise estimations of the 
density of states, and thus it is straightforward to ob- 
tain the microcanonical entropy in computer simulations. 
Since indicative quantities such as transition tempera- 
tures can be quantitatively determined, our method also 
enables experimentally competitive predictions. 

This project has been partially supported by NSF 
DMR-0810223. 
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